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The Landau spectrum of bismuth is complex and includes many 
angle-dependent lines in the extreme quantum limit. The adequacy 
of single-particle theory to describe this spectrum in detail has been 
an open issue. Here, we present a study of angle-resolved Nernst ef- 
fect in bismuth, which maps the angle-resolved Landau spectrum for 
the entire solid angle up to 28 T. The experimental map is in good 
agreement with the results of a theoretical model with parabolic dis- 
persion for holes and an extended Dirac Hamiltonian for electrons. 
The angular dependence of additional lines in the Landau spectrum 
allows to uncover the mystery of their origin. They correspond to 
the lines expected for the hole Landau levels in a secondary crystal 
tilted by 108 degrees, the angle between twinned crystals in bismuth. 
According to our results, the electron reservoirs of the two identical 
tilted crystals have different chemical potentials and carriers across 
the twin boundary have different concentrations. An exceptional 
feature of this junction is that it separates two electron-hole com- 
pensated reservoirs. The link between this edge singularity and the 
states wrapping a three-dimensional electron gas in the quantum 
limit emerges as an outstanding open question. 

Nernst effect | Quantum limit | Landau spectroscopy 

Introduction 

Semi- metallic bismuth is a fascinating solid. Its av- 
erage bulk diamagnetism is larger than in any non- 
superconducting solid, its magnetoresistance is huge and it 
dwarfs any other solid with its Nernst coefficient. During the 
last century, the electronic properties of bismuth were subject 
to numerous studies (See [T||2]|3] for reviews). In particular, 
studies of de Haas- van Alphen effect and Shubnikov- de Haas 
effects led to a precise determination of the components of 
the Fermi surface [U [5], which consists of a hole ellipsoid at 
T-point and three identical electron ellipsoids at L-point[6 . 

Recently, there has been a renewal of interest in the elec- 
tronic spectrum of bismuth at very high magnetic fields in 
the vicinity of the quantum hmitfTl ISl [§l fTOl fTTl . The quan- 
tum limit is attained when the magnetic field is strong enough 
to confine electrons to their lowest Landau level(s). With a 
magnetic field of affordable magnitude, this limit can only 
be attained in low-density metals like bismuth or graphite. 
Emergence of collective effects in a three-dimensional elec- 
tron gas in such a context have been a subject of theoretical 
speculation [12] [13] . 

The measurements of the Nernst response in bismuth up 
to a magnetic field of 12 T led to the observation of giant 
quantum oscillations [14]. The signal was found to sharply 
peak at each intersection of a Landau level and the chemical 
potential. Similar oscillations were also found in graphite [1 5 , 
another bulk semi-metal leading to several theoretical studies 
offering explanations for the occurrence of these giant Nernst 
quantum oscillations 16t (17i il8j . 

The extension of the Nernst data up to 30 T led to the 
observation of additional Nernst peaks at magnetic fields ex- 
ceeding the quantum limit [7 . These peaks were unexpected 
and, given the known structure of the Fermi surface in bis- 
muth, could not be attributed to any Landau level. Moreover, 



angle-resolved torque magnetometry measurements found a 
complex phase diagram for a magnetic field slightly tilted off 
the trigonal axis 8 . These experimental findings raised the 
issue of collective effects in bismuth at high magnetic fields. 
In this context, the high- field phase diagram of bismuth be- 
came subject to a theoretical reinvestigation [T9l I20j . At high 
magnetic fields, the chemical potential continuously shifts to 
preserve charge neutrality. Therefore, quantum oscillations 
are no longer a periodical function of . Further complex- 
ity arises due to an angle-dependent Zeeman energy and cy- 
clotron mass for both electrons and holes 4 . Very recently, a 
fair agreement was found between the results of a theoretical 
model treating the electrons with an extended Dirac Hamil- 
tonian and the experimentally-detected Landau spectrum up 
to 12 T[2T]. The results narrowed down the possible choice of 
theoretical parameters to describe the system in presence of 
moderate magnetic fields but the origin of the additional lines 
remained an open question. 

Several explanations for the origin of additional Nernst 
peaks have been proposed. It was suggested that they had 
been generated by misalignment 20 . But this was ruled out 
by measurements of the Nernst effect in presence of a tilted 
magnetic field TP . It was speculated that surface states could 
give rise to these additional peaks 22 . But both the amplitude 
and the profile of the Nernst peaks point to a bulk origin ^23j. 
Moreover, their complex angle dependence does not follow 
what is expected for a 2D state. Until now, the origin of these 
additional lines on top of the complex Landau spectrum has 
remained a mystery. 

Previous angle-resolved Nernst measurements mapped the 
Landau spectrum in a limited angular window [TTj. Here, we 
present for the first time a map of Landau spectrum extended 
up to 28 T and across the whole solid angle. Comparing the 
experimental map with the results of the theoretical model for 
bulk bismuth allows us to find a solution to this mystery. All 
lines resolved by the experiment match the theoretical expec- 
tations. In particular, the additional lines follow an angular 
variation corresponding to a crystal tilted by 108 degrees as 
expected in case of a twinned bismuth crystal. The results 
suggest that the entire phase diagram of bismuth, at least up 
to 28 T, can be explained within the single-particle theory. 
Surprisingly, however, we find that the agreement between 
theory and experiment can be attained only if one assumes 
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that the chemical potential in the tilted crystal is not set by 
the primary crystal and the two crystals have distinct chemi- 
cal potentials and carrier densities. This implies that carriers 
cannot freely flow across the twinning plane, the interface be- 
tween the two crystals. 

Angle- resolved Nernst effect 

We measured the Nernst effect with a standard one-heater- 
two-thermometer set-up. The Nernst voltage was measured 
with two pairs of electrodes. As the magnetic field rotates, 
we switched from one pair of electrodes to another in order 
to measure the component of electric field in the plane per- 
pendicular to the magnetic field. For each rotating plane, the 
experiment was performed in a dilution refrigerator inserted 
into a superconducting magnet for fields below 12 T and in 
a He^ cryostat inserted in a resistive water-cooled magnet for 
fields between 12 T and 28 T. 

Fig. 1 presents the raw Nernst data as a magnetic field is 
swept in the (trigonal, binary plane). The data for the two 
other planes are shown in the supplementary section. Each 
Nernst peak is associated with the intersection of a Landau 
level and the chemical potential. As seen in the bottom panel 
of the figure, different sets of data obtained at low and high 
magnetic field and with different pairs of electrodes match 
with each other. The peaks associated with the hole pocket 
can be clearly identified and are represented by red solid cir- 
cles in the figure. But other peaks remain, whose origin is 
harder to pin down. 

Theoretical model 

For electrons at L-points, we employed the extended Dirac 
Hamiltonian introduced in ref. i21j. The parameters used 
here is also the same as those in ref. [21], where the fitting is 
carried out below 12 T. 

At high field, however, we have to take into account the in- 
terband coupling between the lowest Landau levels (LLLs) of 
the conduction and valence bands 24 , while it is negligible be- 
low 12 T. At zero field, the conduction and valence bands will 
not be coupled since they have opposite parity. On the other 
hand, at a finite field, the Bloch band picture is no more valid 
due to the complex interband effect of a magnetic field [25", 
which can couple the two LLLs. Thus the coupling should be 
proportional to the field. Such a interband coupling scenario 
actually agrees with the magneto-optical measurements 24 . 

Here we introduce the modified LLL to take into account 
the interband coupling between the two LLLs in the form 




and e{kz) = y^A^ + /i^/cfA/m^, which are obtained by sim- 
plifying the model used in 24 . (/3o = \e\h/meC = 0.1158 
meV/T.) Here, two new parameters are introduced: V ex- 
presses the magnitude of the interband coupling, and V' is 
the correction to the additional g- factor g\ both of which give 
contributions proportional to B. Both V and V' are scalars 
and dimensionless, whereas the longitudinal mass rriz and g' 
are tensors as is given in ref. [21] . 

The separation between the two LLLs decreases with in- 
creasing magnetic field in the low- field region; the separation 
takes its minimum value at Bmin, and then increases in the 
high-field region. We set the parameters as V = 0.15 and 
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V = -0.0625. With these parameters, Bmin ^ 30 T, and the 
minimum value of the separation is about 1 meV (See the sup- 
plementary section for a discussion of the different parameters 
used here and in ref.|24]) . 

For holes at T-point, we employed the same parabolic 
model and parameters used in ref. 21 , except for the value 
of the spin mass along the trigonal direction, Mss = 10000. 
In the previous version of the model, this parameter was 200 
(See TABLE HI in ref. [21]). This minor modification does 
not change the spectrum below 12 T at all, and gives better 
agreement at high fields. 



Crystal structure and twins 

The key point in elucidating the mystery of additional lines 
in the Landau spectrum resides in the rather exotic rhombo- 
hedral crystal structure of bismuth. The unit cell is obtained 
by pulling a cube along its body diagonal. This departure 
from higher symmetry is believed to be triggered by a Peierls 
transition [26]. The distortion can be visualized by considering 
a symmetry- lowering operation on a tetrahedron. In absence 
of distortion, each atom would be at a vertex of a tetrahe- 
dron. At each of these vertices, the angle between two adja- 
cent edges is 60 degrees (Fig. 2a). The tetrahedron remains 
invariant by a 27r/3 rotation around any of the four equiva- 
lent trigonal axes. The structural distortion corresponds to 
a displacement of one of the four vertices, which reduces the 
angle between the three adjacent edges at that particular ver- 
tex to less than 60 degrees. In this distorted structure, out of 
four original trigonal axes, only one survives and the structure 
keeps its threefold rotational symmetry around this trigonal 
symmetry. 

Now, any of the four vertices of the tetrahedron can be 
subject to such a distortion (Fig. 2a). Therefore, bismuth 
crystals can have up to four domains. When at least two of 
such domains are simultaneously present in a crystal, the crys- 
tal is twinned. The atomic positions near a twin boundary in 
bismuth, originally documented by Hall[27 , is illustrated in 
Fig. 2b. A typical bismuth single crystal has many of such 
twin boundaries. Their presence can be used to establish a di- 
rection for trigonal and bisectrix axes and thus determine the 
sign of the tilt angle of the electron ellipsoid 28 . Twin bound- 
aries have been directly detected by STM measurements [29] . 
When a sample is twinned, in addition to the principal crys- 
tal, there is a secondary crystal, which has its trigonal axis 
tilted by 108 degrees respective to the original crystal. The 
two crystals share one binary axis (fig. 2c). 



Landau spectrum 

Fig. 3 compares the theoretical prediction with experimental 
results for three planes of rotation. For each plane of rotation, 
a panel shows both the Nernst peaks resolved by the exper- 
iment and the expected theoretical lines for holes and elec- 
trons. There is no ambiguity for hole peaks. As for electron 
peaks, we have selected Nernst peaks displaying an angular 
dependence close to what is theoretically expected. As seen 
in the figure, the agreement is excellent for holes, and rough 
for electrons. In the case of binary-bisectrix plane (Fig. 4c), 
the experimentally resolved hole peaks are sometimes split be- 
cause of a small residual misalignment. When the field is ori- 
ented perpendicular to the trigonal axis, the Zeeman splitting 
of holes is extremely sensitive to the orientation and suddenly 
vanishes when the field is perpendicular to the trigonal axis 
(See panels a and b of the same figure). 

Given the quality of the agreement between theory and 
experiment in case of holes, it is natural to wonder about 
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the origin and significance of the discrepancy observed in the 
case of electrons. The discrepancy is particularly visible at 
high fields. The theoretical lines for electrons tend to occur in 
magnetic fields significantly lower than the experimental ones. 
This is particularly the case for the angular window around 
the trigonal axis. We failed to find adequate theoretical pa- 
rameters to close this gap between theory and experiment. 



ment between theory and experiment is satisfactory, if one 
assumes that the each of two crystals keeps its own chemical 
potential set by its local environment with respect to the ori- 
entation of the magnetic field. As seen in the lower panels of 
Fig. 5, this means that when the field is near the trigonal axis, 
there is a significant discontinuity both in chemical potential 
and in carrier density across the twin boundary. 



The additional lines 

In Fig. 3, only those peaks which could be attributed either 
to electron or holes were shown. Fig. 4 display additional 
peaks (i.e. those which are clearly not associated with elec- 
trons) together with the hole peaks. The experimental peaks 
are compared with what is theoretically expected for holes 
of a secondary tilted crystal. This secondary crystal shares 
one binary axis with the original crystal, but has its trigonal 
axis tilted 108 degrees off, as depicted in Fig. 2c. As seen 
in the figure, the agreement is good. It is useful to carefully 
examine panel c. As seen in this figure, the additional lines 
merge with the hole lines when the field is aligned along one of 
the three binary. This is the binary axis common to the two 
twinned crystals. When the magnetic field is oriented along 
this binary axis, the spectrum does not show any additional 
features. In this configuration, the two crystals have the same 
crystal axis along the magnetic field and share an identical 
Landau spectrum. 

Previous studies of the high-field phase diagram of bis- 
muth focused on the angular window close to the trigonal 
axis[8l [Til. As seen in Fig. 5, in this range, the results ob- 
tained here are quite similar to those found in a previous study 
of angle-resolved Nernst effect 11 on another bismuth crys- 
tal. In both cases, experiment detects several lines and in 
particular three unexpected high-field peaks when the field is 
oriented along the trigonal axis as originally reported back in 
2007[7j. The additional peaks are not fractional peaks as ini- 
tially assumed[7], but can be indexed as 2^, 1^ and 3^ of a 
secondary tilted crystal as expected in the twinning scenario. 

As seen in Figures 4 and 5, in the case of hole lines, the 
agreement between theory and experiment is impressive, both 
for the primary and the tilted crystal. This agreement has 
been obtained by assuming that the chemical potential contin- 
uously shifts with increasing magnetic field to preserve charge 
neutrality. Now, this field-induced shift in chemical potential 
strongly depends on the relative orientation of the magnetic 
field and the crystal axes [2 1. One can think of two possi- 
bilities in presence of a strong magnetic field: i) the twinned 
crystals share a common chemical potential as they do at zero 
field, ii) The twinned crystals do not share a common chemi- 
cal potential and as the magnetic field is applied, a difference 
between their Fermi energies and their carrier densities grad- 
ually builds up. 

According to our results, the second scenario clearly pre- 
vails. If one assumes that the twinned crystals share a com- 
mon chemical potential set by the primary crystal, then theory 
fails to reproduce the experimentally-resolved spectrum (See 
the supplementary section). On the other hand, the agree- 



Conclusion and open questions 

In a typical bismuth crystal, twin boundaries abound but mi- 
nority domains occupy a small fraction of the total volume of 
the sample [30 . It is therefore quite a surprise that the contri- 
bution of one of these domains is such that the magnitude of 
the Nernst peaks caused by their Landau levels is comparable 
to those of the primary crystal. What sets the magnitude of 
the Nernst quantum oscillations remains an open question. If 
the field-induced wrapping states (the 3D equivalent of the 
edge states of the Quantum Hall Effect [31 ), happen to play 
a role, then a twinned crystal with sufficient depth can gen- 
erate a signal which weighs much larger than the volume it 
occupies. 

The origin and the fine structure of the barrier between 
the crystals with different chemical potential emerges as a 
new subject. Such barriers are common in semiconducting 
heterostructures. But the one under study has a number of 
intriguing peculiarities: the two crystals share the same struc- 
ture and an identical chemical potential in zero magnetic field. 
Moreover, while the carrier density is different across the junc- 
tion, the two systems are compensated and charge neutrality 
is locally preserved. Finally, the discontinuity in the chem- 
ical potential is a consequence of Landau quantification in 
a very anisotropic context. The wrapping states discussed 
by Koshino and co-workers |31j have different Landau indexes 
across the twin boundary. Their length scale is set by the 
magnetic length. On the other hand, the length scale associ- 
ated with charge reorganization at the interface is the Fermi 
wavelength. In our context of investigation, these two length 
scales are comparable in size. 

One conclusion of this study is that the one-particle theory 
can successfully reproduce the experimentally-resolved phase 
diagram of bismuth in a magnetic field as strong as 28 T. 
The fate of the three-dimensional electron gas in bismuth and 
in graphite diverge [32]. In graphite, the system undergoes 
a field-induced thermodynamic phase transit ion [33". This 
phase transition is believed to be a nesting-driven density- 
wave transition 34 , which is expected occur due to the in- 
stability of the quasi-one-dimensional conductor emerging be- 
yond the quantum limit. The absence of a phase transition in 
bismuth raises the following question: what protects bismuth 
from a similar instability? 
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Fig. 1. Top: Nernst signal in a tilted magnetic field in the (trigonal, binary) plane for 
different tilt angles. Bottom: Nernst peaks superposed on a color plot of the same data. 
The angle 6 represents the angle between the magnetic field and the trigonal axis. Those Nernst 
peaks, which are clearly associated with holes are represented by solid red circles and others 
by white triangles. The high- (low-)field data was obtained at 1.5 (0.49) K. Digits represent 
orbital indexes for hole Landau levels. They become spin-split (zb) when ^ 90°. 
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Fig. 2. a) Structural distortion in bismuth: The crystal structure of 
bismuth can be assimilated to a distorted tetrahedron. The tetrahedron can be distorted by 
pulling along any of the four different trigonal axes. At the corresponding vertex, the angle 
between each pair of adjacent edges becomes lower than 60 degrees, b) Sketch of a twin 
boundary in bismuth after [27]: Open circles represent atoms in the undistorted 
crystal. Solid circles represents atomic positions in the twinned crystal. The angle between the 
trigonal axes of the two crystals is about 108 degrees. The common binary axis for the two 
crystals is perpendicular to the plane of the figure, c) Conventions used in this 
text: The green plane and axis are the original trigonal plane and axis. The blue ones are 
tilted planes and axes. The two structures share one binary axis. 
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Fig. 3. The Landau spectrum for holes and electrons according to theory (solid lines) and 
experiment (symbols). The three panels correspond to the three rotating planes. Peaks resolved 
by experiment and identified as arising from holes (red open circles) or electrons (green open 
triangles) are displayed. 
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Fig. 4. The Landau spectrum for holes in the original and tilted crystal according to theory 
(solid lines) and experiment (symbols). The three panels correspond to the three rotating 
planes. Additional peaks resolved by the experiment are represented by blue open squares. 
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Fig. 5. Top: Experimental (a) and theoretical (b) Landau spectrum for field oriented close 
to the trigonal axis. The lines are indexed according to their orbital and spin quantum number. 
The additional lines correspond to tilted hole lines. Bottom: The evolution of the carrier density 
(c) and chemical potential (d) for the primary (Dq) and tilted(Dt) crystals when the magnetic 
field is parallel to the trigonal axis. 
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PACS numbers: 
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I. EXPERIMENTAL DETAILS 



We measured the Nernst effect with a standard one-heater-two-thermometer set-up, which 
was rotated in a magnetic field using Attocube piezoelectric rotators. All three perpendicular 
rotating planes were scanned in order to extract the map for the whole solid angle. In the 
main text, the raw data were presented for one plane of rotation. We present in Fig. 1 and 
Fig. 2, the data for two other planes of rotation. 

II. THEORETICAL DETAILS 

The model used here is slightly different from the one used in ref.[l]. The magnetic field 
dependence of the lowest Landau levels (LLLs) for the present parameters {V — 0.15 and 

— —0.0625) are shown for B \\ binary in Fig. 3 (a). The origin of the energy is at the 
center of the band gap, while in Fig. 3 of ref. [1], it is taken at the bottom of the conduction 
band. 

The current model takes into account an interband effect induced by magnetic field, which 
can couple the two Lowest Landau levels (LLLs). In the present case, the separation between 
the two LLLs takes its smallest value at ^min ~ 30 T, and the LLL in the conduction band 
does not cross the Fermi energy up to 45 T. 

The result with V — 1.26 and — 0.253, which are derived from the parameters used 
in ref. [2], is shown in Fig. 3 (b). The result agrees well with Fig. 2 of ref. [2], even though 
our model is slightly different from theirs. (Note that the result shown in Fig. 2 of ref. [2] 
is not the result obtained experimentally, but is the result obtained theoretically.) In their 
case, the separation takes its smallest value at ^min ~ 10 T, and rapidly increase above 10 
T, due to the large values of V and V\ 

As seen in Fig. 3 (b), the LLL in the conduction band crosses the Fermi energy at 
around 20 T, and the electron pocket becomes completely empty above this field. This 
crossing would have generated an additional peak. However, such a peak has not been 
observed in our Nernst experiments up to 28 T. Therefore, the current experimental results 
favor the situation sketched in Fig. 3 (a). 
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FIG. 1: Supplement: Top: Nernst signal in a tilted magnetic field in the (trigonal, bisectrix) plane 
for different tilt angles. Bottom: Nernst peaks superposed on a color plot of the same data. 
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FIG. 2: Supplement: Top: Nernst signal in a tilted magnetic field in the (binary, bisectrix) plane 
for different tilt angles. Bottom: Nernst peaks superposed on a color plot of the same data. 
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FIG. 3: Supplement: Magnetic field {B \\ binary) dependence of the lowest Landau levels (Oe-) 
and the Fermi energy {Ep) for (a) V = 0.15 and = —0.0625 (present work), and (b) V = 1.26 
and = 0.253 (ref. [2]). The dashed lines are the results without considering the interband 
coupling between the two lowest Landau levels. 

III. FAILURE OF THE COMMON CHEMICAL POTENTIAL SCENARIO 

In the manuscript, we have compared theory and experiment in the case of diflFerent 
chemical potentials for the primary and tilted crystals (Fig. 4 and Fig. 5 in the main 
manuscript). In Fig. S4, we show the results of theoretical calculations for one rotating 
plane assuming a common chemical potential for both crystals. As seen in the figure, there 
is a significant discrepancy between theory and experiment. The failure of theory by a very 
large margin in this scenario compared to its success when the chemical potential assumed 
to be diflFerent backs one of our main conclusions. 
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FIG. 4: Supplement: The Landau spectrum for holes in the original and tilted crystal according 
to theory (solid lines) and experiment (symbols) for the (trogonal, bisectrix in case of a single 
chemical potential. Additional peaks resolved by the experiment are shown by blue empty squares. 
They do not match the theoretical lines. 
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